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The small-scale turbulent dynamo in the high Prandtl number regime is 



0^ ' described in terms of the one-point Fourier space correlators. The second or- 

der correlator of this kind is the energy spectrum and it has been previously 

> 

' studied in detail. We examine the higher order A:-space correlators which 

' contain important information about the phases of the magnetic wavepackets 

o 
cn 
o 



and about the dominant structures of the magnetic turbulence which cause 
intermittency. In particular, the fourth-order correlators contain information 
about the mean-square phase difference between any two components of the 
magnetic field in a plane transverse to the wavevector. This can be viewed as a 
measure of the magnetic field's polarization. Examining this new quantity, the 



^ ' magnetic field is shown to become plane polarized in the Kazantsev-Kraichnan 



model at large time, corresponding to a strong deviation from Gaussianity. 
We derive a closed equation for the generating function of the Fourier cor- 
relators and find the large-time asymptotic solutions of these correlators at 
all orders. The time scaling of these solutions implies the magnetic field has 
log-normal statistics, whereas the wavenumber scaling indicates that the field 

is dominated by intermittent fluctuations at high k. 
47.27.-i, 95.30.Qd 
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I. INTRODUCTION 



In some astrophysical applications such as the interstellar medium and protogalactic 
plasmas the kinematic viscosity is greater than the magnetic field diffusivity by the factor 
10^"^ to 10^^ [1)2]. In these situations, there is a vast scaling range where the magnetic field 
has a smaller characteristic length-scale than the velocity field. The dynamo process of the 
stochastic stretching and amplification of the magnetic field can be studied in this case in 
terms of the statistics of Lagrangian deformations [3,4]. It can be pictured as a collection 
of magnetic field wavepackets, each moving along a fiuid particle path and being distorted 
by the local strain.^ Such a regime of smooth velocity fields is similar to the Batchelor's 
regime found in the related problem of passive scalar advection. The Batchelor's regime has 
a long history of study, some of the more recent advances can be found, for example, in [2,4] 
(and references therein). A common further simplification of the problem is to assume that 
the local strain matrix is a Gaussian white noise process, this is commonly known as the 
Kazantsev-Kraichnan model [5,6]. However, this assumption is not always necessary, and 
some results have been shown to be universal for a broader class of stochastic fiows [3]. 

The first analysis of the dynamo problem were based in Fourier space [5-7] , with par- 
ticular emphasis being placed on the second order moment, corresponding to the energy 
spectrum. It was realized, however, that the Fourier space and coordinate space equations 
have a similar structure. Recent studies of the turbulence intermittency (in both the passive 
scalar and dynamo problems) have focussed on coordinate space moments of second order or 
higher. This approach was motivated by a feeling that this description is more natural and 
can give more information than the /c-space moments [3,4,8]. As a result, the only /c-space 
correlator that has been seriously studied to date is the energy spectrum and there is no 
theory describing the higher Fourier space correlators. 

In the present paper, we turn our attention back to Fourier space and consider the 



The second and higher derivatives of the velocity field can be ignored in this case. 
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one-point moments of the magnetic field Fourier transform. Although the problem under 
investigation here is the small-scale turbulent dynamo, we would like to give several rea- 
sons why the description of the Fourier space moments is important and why it should be 
developed in a broader turbulence context. 

(i) The presence of singular structures in turbulence is known to affect the scaling of the 
structure functions. The most well-known examples of this being the /3 and multi-fractal 
models [9]. However, some coherent structures are singular in Fourier space rather than in 
coordinate space, and therefore can be detected by an investigation of the Fourier moments. 
A simple example of structures that are singular in Fourier space (and regular in coordinate 
space) is a sea of vortex filaments in 2D turbulence. The layered pattern of the vortices in 
coordinate space corresponds to a ID curve in Fourier space. ^ 

(ii) To date, the only Fourier moment that was studied in detail is the second order 
moment that describes the distribution of the turbulent energy over wavenumbers. However, 
there exist other Fourier moments that have a clear physical meaning and describe important 
properties of the turbulence. One of these objects is related to the fourth order moments 
and has a physical meaning of the mean turbulence polarization. It will be introduced and 
studied in this paper. 

(iii) Phases of the Fourier modes can be dealt with directly, and therefore the validity of 
the random phase assumption can be examined. 

(iv) Finally, in some cases a Fourier space analysis is the only way to have a treat- 
able problem due to the introduction of non-localities, for example via a pressure or wave 
dispersion. Note that the dynamo and passive scalar systems do not fit into this class of 
problems. However, the methods developed in the present paper will be applied to the 



^The Fourier transform should be taken over a local box, which is large enough to fit many layers 
of filaments, but is smaller than the large-scale vortices in this case. This will be discussed in the 
next section. 
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non-local Navier-Stokes equation^ (which involves a pressure term) in our next paper [10]. 

II. ONE-POINT AND TWO-POINT CORRELATORS OF FOURIER 

AMPLITUDES 

To avoid any confusion, we should state that the Fourier transforms in this paper are 
performed over a finite box of a size much greater than typical length-scale of the magnetic 
turbulence, but much less than the scale of the advecting velocity field. It is only for finite 
box Fourier transforms that one-point correlators are well defined objects. Secondly, the 
center of each finite box is not fixed. The box center coordinate dependence of the Fourier 
transforms is a useful measure of the slow variability of the large-scale magnetic turbulence. 
In this paper we make each box move in unison with a fiuid particle, and thus any coordinate 
dependence is replaced by a time dependence along a given fiuid trajectory. 

One should be aware that in this case the finite box size does not produce a discrete 
spectrum as the magnetic field is in general not periodic within the box. On the other hand, 
the lack of periodicity can produce unphysical power-law tails in the Fourier transforms at 
(large) wavenumbers that are greater than one over the box length. Slow decay of these tails 
can be undesirable in numerical simulations and other applications. However, this problem 
can be easily cured by using a smooth filter instead of the box, for example by using Gabor 
transforms instead [11]. In this presentation the box (or filter) shape is actually not an issue. 
We will present the results in terms of finite box Fourier transforms. However, performing 
the analysis using Gabor transforms will in fact produce identical results. 

Let us consider the general two-point correlator of the Fourier transformed magnetic field 
components 



^Such an equation is the basis of Rapid Distortion theory which will be developed for the case of 
a stochastic strain. 
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where 11,12, I'^n and ji, J2, ■ ■ ■ ,jn take the values 1, 2 or 3 corresponding to the three 
components of the magnetic field, and n is an arbitrary natural number. Let us assume 
the turbulence is quasi-homogeneous. That is, the correlator (i?i(xi)i?j(x2)) (and similar) 
depend only on Xi + X2 as slowly as the advecting strain field. On the other hand, such 
correlators decay in Xi — X2 at distances much less than the box size L. In this case 

(5,,(k05,,(ki) . . .5,„(k05,,(k2)5,,(k2) • • •%(k2)) = (2) 
{B,,{ki)B,,{ki) . • •S,„(ki)%(-ki)S,,(-kO • • • %(-ki)) /(ki + k2), 

where /(k) = -^Ilf^^ ^'"^^^^''' is the Fourier transform of the filter function and d is the number 
of space dimensions. Here we see that any two-point correlator is fully determined by the 
one-point correlators 

^n..,.-..„;.i,..,-....(k) = (i?n(k)5..(k) • • •S,„(k)i?,,(-k)i?,,(-k) • • •%(-k)). (3) 

This simple observation means we can avoid lengthy derivations dealing with two-point ob- 
jects directly, and instead can obtain results from the one-point correlators that are them- 
selves far easier to deal with. Turbulence isotropy and the divergence-free condition further 
narrow the class of possible one-point (and two-point) correlator tensors we need to consider 
(see for example [12]) 

4,i2,-,in;jlJ2,-,Jn(k) = Yl ^''sik)Dl,,l,Dl,^U " " " A„_l,i„, (4) 

pairs 

where 

A. - - (5) 

and the summation is over the set of all possible pairings (/i, I2), {I3, U), - ■ ■ , {In-i, In) of the 
indices i and j from the left hand-side of (4). For each pairing, the index s is equal to the 
number of pairs that consist only of i's, ^ for example (^1,^3) or {i%,i2)- One can see that 



^There are obviously going to be the same number of pairs that consist only of j's. 
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any correlator of order 2n can be expressed in terms of N — Max(s) + 1 = Int(n/2) + 1 inde- 
pendent functions ^'^{k). There is only one such function for the second order correlators, 
two functions for the fourth and sixth orders, and three functions at the eighth and tenth 
orders. 

Instead of we can choose another set of independent functions, in particular the 

following set of correlators 

= (|B(k)|(2"-^*)|B(k)Y')- (6) 

These correlators are linear combinations of $"(A;)'s, for example at the fourth order we 
have 

=2(d-l)$[,'^ + (d-l)2$f). 

Below, we will derive a technique that allows us to deal with the fundamental set of corre- 
lators (6). 

III. THE MODEL 

Let us start with the equation for the magnetic field B(x, t), 

atB + (u-V)B = (B-V)u + fi:AB, (7) 

where u is the velocity field, which we assume has a much slower spatial variation than the 
magnetic field B, and k is the magnetic diffusivity (determined by the fiuid conductivity). 
We will make use of the previously discussed finite box Fourier transforms. Each box has 
sides of length L which is chosen to lie between the length-scales associated with B and u, 
namely Lb and L„. The box Fourier transform is defined as 

B(k,x,i) f B(r, e'''- dr, (8) 

J box 
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where x is the coordinate of the box center. Applying this Fourier transform to (7) we have 
(with accuracy up to the first order of the scale separation parameter e = Lb/L^ <^ 1) 

dtBm + UiViBm = (Tijh djBm + (JmiBi - Kk^B^, (9) 

where aij — V jUi is the strain matrix and the operators Vj and di correspond to derivatives 
with respect to Xi and ki respectively. In considering a fluid path determined by x(t) = u, 
and thus noting that aij{x.,t) aij{x.{t),t) and B(k, x, B(k, x(i),i), this equation 
becomes^ 

dtBm = (Jijki djBm + (JmiBi - Kk^Bm- (10) 

It should also be noted that the strain a^j, taken along a fluid path, only enters this equation 
as a given function of time. To complete the model one has to specify this dependence or 
postulate the strain statistics. In what follows we choose the strain matrix to be Gaussian 
such that 

= n (^Aij - ^Sij^ , (11) 

where Ay is a matrix with statistically independent elements that are white in time 

{Aij{t)Aki{0)) = SijdkiS{t). (12) 

This choice of strain matrix ensures incompressibility and statistical isotropy. Indeed, 

{aij{t)aki{0)) = n (^dik6ji - ^^ijdki^ d{t). (13) 

However, one must be aware that this is not the only way to satisfy the incompressibility 
and isotropy conditions; for example [3,4,13] have chosen 

{(Tij{t)aki{0)) ^n[{d+ l)SikSji - SilSjk - SijSki] S{t). (14) 



^Hereafter, we will the drop hat notation on B because only Fourier components will be considered. 
Also, we will not mention explicitly dependence on the fluid path and simply write B = B(k, t). 
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More generally, there is an infinite one-parametric family of possible strains. The strain 
matrices that are Gaussian and white in time correspond to the Kazantsev-Kraichnan model, 
a natural starting point for an analytical analysis because of its simplicity. However, it should 
be noted that some results remain universal in the case of other smooth velocity fields found 
in a much wider class of statistical models [3,4,13]. In our future work we will investigate 
the behaviour of the Fourier space correlators in the case of more general statistics. 



IV. GENERATING FUNCTION 

Let us consider the following generating function 

Z{\ a, /3, k) = (e^|B(k)P+aB^(k)+/3B^(k)^^ (15) 

where the overline denotes the complex conjugation. This function allows one to obtain any 
of the fundamental one-point correlators (6) via differentiation with respect to A, a and (5 

^;^= (|B(k)|(2"-4-)|B2(k)|2^) = bf^^^'^^S^a^zl . (16) 

L J A=a=/3=0 

Differentiating (15) with respect to time and using the dynamical equation (10), we have 

Z = kidjiaijE) + \{(Tmi{BmBi + BiB^)E) + 2a{amiB^BiE) + 2p{amiB^BiE) (17) 
-2ne{[X\B{k)\^ + aB^{k) + f3B\k)]E), 

where 

^ ^ gA|B(k)|2+aB2(k)+/3B'(k)_ (^g) 

To find the correlators on the right hand-side of (17), one needs to make use of the Gaus- 
sianity of the strain matrix aij and perform a Gaussian integration by parts. We then use 
the whiteness of the strain field to find the response function (functional derivative of Bi 
with respect to (jy). Finally, one can use the statistical isotropy of the strain, so that the 
final equation involves only A; = |k| and no angular dependence of the wavevector. This 
derivation is discussed in more detail in the Appendix. Here, we just write the final result 



(19) 



- 2Kk'^VZ, 



where the k, a, (5 and A subscripts in Z denote differentiation with respect to k, a, (5 and A 

respectively, and 



V = \dx + ada + Pd)3. 



(20) 



Actually, one can reduce the number of independent variables in this equation by taking into 
account that Z depends on a and /3 only in the combination r] — aj3. The easiest way to see 
this is to consider a Taylor series expansion of Z from (15); any term that contains a and (5 
in a different combination will be zero because of the quasi-homogeneity of the turbulence. 
Thus, we can write 



1 - ^ k^Zkk + i(-4P + d^- l)kZk 



(21) 



+ {2d - 6)VZ + 4[1-1]V'Z + 2{X^ - 47j){Z, + tjZ,, - Z^^) 



2Kk^VZ, 



where 



V = Xdx + 27]d^. 



(22) 



In what follows, we will restrict our consideration to the 3D case {d — 3) in which case the 
equation for Z reduces to 



Z = ^ [eZkk + (4 - 2V)kZk 

+AV^Z + 3(A' - Ar)){Zr, + r]Zr„ - Zxx)] - 2Kk''VZ. 



(23) 



V. ENERGY SPECTRUM 



Let us now consider the energy spectrum of the magnetic turbulence, given by the second 
order correlator 
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E{k,t)^^l^{\B{k)\')^[d,Z],^^^,. (24) 

Differentiating (23) with respect to A and taking the result at A = 7y = we have 

E^^{k'^Ekk + 2kEk + 4:E)-2nk'^E. (25) 

This equation for the evolution of the energy spectrum was first obtained by Kazantsev [5] 
and Kraichnan and Nagarajan [6]. Kazantsev analyzed an eigenvalue problem associated 
with this equation which allowed him to obtain the growth exponents of the total magnetic 
energy. Numerically, the energy spectrum was studied by Kulsrud and Anderson [1] who 
gave a detailed description of the /c-space evolution of this spectrum. Recently, Schekochihin, 
Boldyrev and Kulsrud [2] presented the complete solution of (25) obtained by the use of 
Kontorovich-Lebedev transforms (KLT). Note that this integral transform approach has an 
advantage over the Kazantsev's eigenvalue analysis in that it allows us to obtain not only the 
growth exponents, but also the power-law prefactors, of the large time asymptotic solutions. 

As they will be of use later in this presentation, let us briefly review the previous results 
for the energy spectrum before we discuss the higher order moments. Using the following 
substitution [2] 

E = e^^*/^ k-'/' cPik/ka, , t) k,„ = (26) 

one can reduce (25) to 

3 • 

-0(p, t) = p^(ppp + p0p - pV, (27) 

where p — k/kd^. At scales much greater than the dissipative one, p <S 1, there is a perfect 
conductor regime for time t <C (In q) ^ (where g <C 1 is the mean wavenumber of the initial 
condition). Thus, in this regime the last term in (25) can be neglected. By changing to 
logarithmic coordinates and a moving frame of reference one can transform this equation 
into a heat equation. For 1/VL the solution of which is just the Green's function, which 
gives [1,2,14,15] 

= const r^'^ e-^^, (28) 
10 



where the constant is fixed by the initial condition. This solution describes a spectrum with 
an expanding k~^/'^ scaling range. At t ~ (Ing)^ the front of this scaling range reaches the 
dissipative scales. To solve (25) in this case, we note that the right hand-side of this equation 
is just the modified Bessel operator and by using KLT one immediately obtains [2] 

const / dsssmh{'Ks)Kis{p)Kis{q)e * *, (29) 
Jo 

where Kis is a MacDonald function of imaginary order. Again this solution is given by the 
Green's function only because the condition that i » 1/Q is obviously satisfied if g <S 1. 
For time t ^ (Ing)^, the function e~^^^ is strongly peaked at t = and integration of (29) 
gives 

= const t-^/^Ko{p). (30) 

One should note that although the Ko{p) shape is somehow predicted by the Kazantsev's 
eigenmode analysis [2,5,14,15], the t^^'^ factor can only be obtained by solving the full initial 
value problem. For p -C 1, Kq{p) ^ — Inp, which means that at large times the scales far 
larger than the dissipative one are affected by diffusion via a logarithmic correction 

E{k) ~ k-^/^ Hkdjk). (31) 

The energy evolution gives an important, but incomplete picture, of the dynamo process. In 
particular, it does not capture the existence of small-scale intermittency and does not allow 
us to predict the type of coherent structures dominating the turbulence at large time. To deal 
with these issues one has to study higher order correlators. Higher one-point correlators in 
coordinate space were studied in [3] and used to predict the shape of the dominant structures. 
Below, we proceed to study the higher fc-space correlators. In particular, this will lead to 
the discovery of a new quantity of interest, corresponding to the mean polarization of the 
magnetic turbulence. 
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VI. 4TH-ORDER CORRELATORS, TURBULENCE POLARIZATION AND 

FLATNESS 



There are two independent 4th order correlators 

S{k,t) = v^(^) = (|B(k)r) = [Z,,],^^^, , (32a) 

T{k,t) = = (|B2(k)n = [Z,],^^^, . (32b) 

Differentiating (23) twice with respect to A and taking the result at A = = we have 

S^^ (eSkk + lOS + 6T) - Ank'^S. (33) 

Now, differentiating (23) with respect to rj and taking the result at A = ?7 = we get 

T = - {k'^Tkk + 4r + 125) - Ank'^T. (34) 
3 

Equations (33) and (34) make up a complete system for 5" and T and can be solved exactly 
in the general case. Observe that there is a closed equation for = 5 — T 

W = - {k'^Wkk - 2W) - Ank'^W. (35) 
3 

Before solving this equation, let us examine the physical meaning of W by writing it as 

3 3 

W = 4^([S(i?,5,)]^) = 4^(|i?,f |i?,|2sin^(0, - 00) > 0, (36) 

where $j denotes the imaginary part and and 0; are the phases of the components Bj and 
Bi respectively. We see therefore that W contains information not only about the amplitudes 
but also about the phases of the Fourier modes. In particular, W = corresponds to the 
case where all Fourier components of the magnetic field are plane polarization, li W ^ 
then other polarizations (circular, elliptic) are present. This is the case for example for a 
Gaussian field where one finds W — > 0. On the other hand, the smallness of the 

phase differences in W can be overpowered by large amplitudes. Therefore, a better measure 
of the mean polarization would be a normalized W, for example 
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P = W/S. (37) 

Defined in this way, the mean turbulence polarization is an example of an important phys- 
ical quantity which can be obtained from the one-point Fourier correlators, and that is 
unavailable from the (one-point or two-point) coordinate space correlators. 

In the perfect conductor regime, when the diffusivity term in the equations (33) - (35) 
can be ignored, the solution for W can be obtained in a similar manner to the previous 
energy spectrum analysis. Namely, by reducing (35) into a heat equation via a logarithmic 
change of variables and passing into a moving frame of reference. The solution therefore is 

W^Wo r V^e-^^"*/^ k'/'e-"^ , (38) 

where Wo is a constant which can be found from the initial conditions. We see that W 
develops a k^^'^ scaling range which is cut-off at low and high k by exponentially propagating 
fronts. Within this scaling range, W decays exponentially in time. 

Given W, one can also find S by representing it as 5" = V+cW and choosing the constant 
c such that the equation for V is closed. This gives c — 3/7 and 

S = r V2 k^'e-"^ (Vo e^^^*/^ + ^Wo e"^"*/^^ , (39) 

where Vq is another constant which can be found from the initial conditions. For t ^ 

the second term in the parenthesis should be neglected, and we have the following solution 

for the mean turbulence polarization 

P^W/S^ (40) 

In the perfect conductor regime we see the mean polarization tends to a value that is 
independent of k and decays exponentially in time. This means that all Fourier modes of 
the magnetic field eventually become plane polarized. Recall that such turbulence is very far 
from being Gaussian, remembering that the mean polarization of a Gaussian field remains 
finite.^ 



^That is, elliptic and circular polarized modes are present. 
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It is also easy to obtain solutions to W for the diffusive regime. Indeed, following the 
example of the energy spectrum we make use of the substitutions 

M^ = e-^"*/^fcV2^(fc/A;,J, (41) 

y ^g2im/4^i/2^^^/^^j^ (42) 

to transform the governing equations for W and V into a form similar to (27). In each case 
we can solve the equation for using KLT and find that for times t » (Ing)^ 

W ~ Wok^/H-^/^e-^'''/^Ko{k/kd,), (43) 

V ~ Vok^'H-"^e'''^''^Ko{k/kd,). (44) 



Correspondingly, 5* has the solution 



S ~ k^'H-^"^ 



Mk/kd,), (45) 



while importantly we find that the normalised polarisation P = W/ S behaves identically in 
both the diffusive and perfect conductor regimes. 

Therefore, in the diffusive regime W continues to decrease in time with the same ex- 
ponential rate as perfect conductor case and has the same prefactor as the energy 
spectrum. Thus, by the time the diffusive regime is achieved W can be essentially put 
equal to zero. The fact that the polarization becomes plane has quite a simple physical 
explanation. Indeed, a magnetic field wavepacket of arbitrary polarization will be strongly 
distorted by stretching. The stretching being strongest along the direction of the dominant 
eigenvector of the Lagrangian deformation matrix (corresponding to the greatest Lyapunov 
exponent). Such a stretching will make any initial "spiral" structure fiat at large time, with 
the dominant field component lying in a plane passing through the eigenvector stretching 
and wavevector directions.^ 



^Of course, the incompressibility of B(k) ensures that it is perpendicular to k. 
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Another measure of intermittency in turbulence is the flatness which can be defined in 
/c-space as F — S/E'^. For large time, in the perfect conductor regime 

F~iV2eiint/4^3/2_ (4g) 

We see that the flatness grows both in time and in k which indicates the presence of small- 
scale intermittency. Such an intermittency can be attributed to the presence of coherent 
structures in A;-space. In the diffusive regime, in contrast to the polarisation the behaviour 

of the flatness is modifled and at large time one finds 

P ~ 7,3/2-3/2 iim/4 Ko{k/kd,) , . 

For small k, we again have a region of /c^/^ scaling but now with a log correction arising from 
the MacDonald functions. For large k, below the spectral cut-off, the additional MacDonald 
functions act to heighten the fiatness. That is, the introduction of a finite diffusivity actually 
increases the small-scale intermittency. In the next section we will investigate the correlators 
of all orders. 



VII. LARGE-TIME BEHAVIOR OF HIGHER CORRELATORS 

The observation in the last section, that there is a dominant field component, allows us 
to predict that for large time |B|^ |B^p in each reahzation, that is Zx\ ~ Z^/s. Therefore, 
the property Z\\ = Zap, if valid initially, should be preserved by the equation for Z. Indeed, 
let us differentiate the equation (23) twice with respect to A and subtract it from the same 
equation differentiated with respect to a and /3. This gives the following closed equation for 
the combination w — Zx\ — Zap 

[k'^Wkk - 2kVwk + 4X>^w (48) 
+16{Vw + w) + 3{dxx - dafsHiX"^ - 4a/?)w}] - 2Kk'^{Vw + 2w). 

We see that if w = at i = 0, then it will remain identically zero for all time. Thus, we can 
consider a class of solutions of (23), corresponding to large-time asymptotics of the general 
solution, such that Zx\ — Zap- Assuming this equality in (23) and putting 77 = we have 
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Z = -{k^ZkH + (-2AaA + 4)fcZfc + ^{\dxfZ\ - ^Kk^XdxZ. (49) 

o 

Let us consider a solution to this equation that is formally represented as a series in A (for 
example via a Taylor series): 

Z=l + J]^*(-)(A;,i). (50) 

n=l 

Here, the function \1'*^") is the correlator of order 2n. Wc omitted the lower index in ^'^"^ 
because the correlators corresponding to different lower indices are identical in this case. 
Substituting (50) into (49) we have the following equation for these correlators 

^ %kH^^l + (-2n + + 4n'*(")] - 2«A;'n*("). (51) 

Note that forn = 1 this equation agrees with the Kazantsev equation for the energy spectrum 
(25). Moreover, by the substitution 

^(n) ^ g(n-l/2)(n+3/2)m ^n-3/2 ^{kjk^^^t), (52) 

one can reduce (51) to an equation for 0, similar to (27), that is independent of n. We can 
therefore immediately write down the general solution for the correlators at any order n. In 
particular, in the perfect conductor regime l/f2 <C t <C (Ing)^ we have 

- COnst(n) g{n-l/2)(n+3/2)m ^n-3/2^ (53) 

and in the diffusive regime t 3> (Ing)^ 

= const(n) r=^/2e("-V2)("+3/2)nt ^n-3/2 Ko{k/kdJ. (54) 

We see that the main effect of the dissipation on all the moments (including the energy 
spectrum) is the prefactor changes from t~^/^ — t~^/^ (but without a change in the expo- 
nential growth rate), further the KQ{k/kd„) form-factor corresponds to a log-correction at 
k <ti I and an exponential cut-off at larger k.^ It is this exponential cut-off that causes 



^Note that the higher the order of the correlator the earlier the spectral cut-off. 
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the exponential growth of the mean magnetic energy [5] and higher x-space moments of the 
magnetic field [3] to change. Indeed, simply integrating ^'^"^ over k with a cut-off at k = kd„ 
(and ignoring the pref actor and log corrections) one recovers the Kazantsev growth rate of 
the magnetic energy in the dissipative regime [5] . Such an explanation was previously given 
in [1]. 

VIII. PHYSICAL INTERPRETATION OF THE SCALINGS 

Let us analyze the physical origins of the scalings obtained in the previous section because 
this can give us an indication to whether the same scalings should be expected for more 
general strain statistics. 

Let us consider the expression (26) for the n-th order correlator and re-write it in the 
form C{n) ■ exp{{2n'^ + 2n)Vtt)k'^ ■ P{k, t), where P(/c, t) is a universal function. The physical 
meaning of various terms in this expression can be easily analyzed. It follows from the 
Central Limit theorem that the large time statistics of the Lyapunov exponents is Gaussian 
with a dispersion D{t) ~ y/i. It follows from time reversal invariance that the average 
values of the Lyapunov exponents are A, 0, —A. Hence, in the large time limit, B{t) ~ 
exp{Xt), where p.d.f.{t) ~ exp{— ^^M-^ )• Therefore, the statistics of the magnetic field is 
log- normal and (5^") ~ exp{{2\n^ -\- 2An)t). For the statistics chosen in the present paper 
A = A = fi. Substituting these values into the last expression for {B"^"-), we restore the 
correct n-dependence of the exponential growth rate of Note that terms of order in 
the growth rate are due to the Gaussian nature of the fluctuations of A around A, while the 
terms of order n are due to the fact that A 7^ 0. It is interesting that the magnetic field 
has log-normal statistics in both the perfect conductor and the dissipative regimes. For the 
coordinate space moments of B, persistence of the log- normality was emphasized in [16], 
although the e^"-^^ dependence was established earlier in [3]. An equivalent result for the 
random matrices can be traced back to [17] (see also [4]). 

The /c-dependence of the magnetic field correlators is also very important because it gives 
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us information about the dominant structures in wavenumber space. Suppose that initially 
the magnetic turbulence is isotropic and concentrated in a ball centered at the origin in 
wavenumber space. For each realization such a ball will stretch into an ellipsoid with one 
large, one short and one neutral dimension. One can visualize this ellipsoid as an elongated 
flat cactus leaf with thorns ahgned to the direction of the magnetic fleld. Note that in this 
picture one component of the magnetic field (transverse to the cactus leaf) is dominant, 
this is captured by the fact that the polarization W introduced in this paper tends to zero 
at large time. Another consequence of this picture is that the wavenumber space will be 
covered by the ellipsoids more and more sparsely at large k, this implies large intermittent 
fluctuations of the magnetic fleld Fourier transform. These fluctuations can be quantifled by 
the flatness F which was shown in (46) to grow as /c^/^, a clear indication of the small-scale 
intermittency. 

To investigate the Kraichnan-Kazantsev model based dynamo problem further, a set of 
numerical experiments have been performed to, among other things, investigate the sen- 
sitivity of these analytical results to changes in the strain statistics. The details of this 
investigation can be found in [14,18]. To help visuahse the cactus leaf structures described 
above we will briefly include some snapshots of the conflguration of a wavepacket ensemble 
in wavenumber space. As above each wavepacket is initially randomly distributed on a unit 
sphere in k-space with a randomly orientated magnetic fleld that lies in a plane perpendicu- 
lar to the wavepacket's wavevector, tangent to the sphere. Figure 1 shows the real magnetic 
fields of a set of 500 wavcpackcts that have been subjected to two different reahsations of the 
strain matrix.^ In the left-hand figure it is clear the magnetic field in this realisation is far 
from being plane polarised, the magnetic vectors still being predominantly random in their 
orientation at this given point in time. In contrast, in the right-hand figure, which has also 
been taken at the same point in time, we see the ellipsoid has become very elongated and 



The imaginary magnetic fields are qualitatively similar and have therefore not been included. 
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the magnetic field appears plane polarised. In contrast, figure 2 shows the same two strain 
realisations but with a finite diffusivity. The right-hand figure in particular demonstrates 
why, in the diffusive regime, an ellipsoid will cover /c-space more sparsely due to the decay 
of its magnetic field at its tips. This is the reason why spectral flatness increases in the 
diffusive regime. 

IX. CONCLUSION 

In this paper we introduced a description of the small-scale magnetic turbulence in terms 
of the one-point correlators of the Fourier amphtudes. Prom the classic works of [5-7] to 
more recent analysis (see for example [2] and references therein) , the only correlator of this 
kind that was studied is the energy spectrum. The higher order correlators have only been 
considered in coordinate space [3,4]. In this paper we have considered the Fourier correlators 
of all orders and showed that they contain important additional information about the 
turbulence, which is unavailable from a similar analysis of one-point or two-point coordinate 
space correlators. In particular, the fourth order Fourier correlators carry information about 
the mean polarization of the magnetic field modes. We showed that this polarization becomes 
plane for the Kraichnan-Kazantsev model. The /c-scahng of the higher correlators allows us 
to determine the structures in Fourier space responsible for the intermittency, which for the 
Kraichnan-Kazantsev dynamo turns out to be elongated ellipsoids centered at the origin. 
The time scahng of the higher correlators allows one to conclude that the magnetic field 
has a log-normal statistics, although the same information is contained, and was established 
before, in analysis of the coordinate space correlators [3,4,16]. 

Finally, we would like to discuss an interesting connection between this work and the 
recent work of Schekochihin et al. [8,16]. In particular the connection between the statistics 
of the magnetic field curvature studied by Schekochihin et al. and the magnetic polarization 
measure introduced in this present paper. Schekochihin et al. found that that the curvature 
of the magnetic field decreases, corresponding to folded and strongly stretched structures. 
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This agrees with our results that the Fourier modes of the magnetic field tend to a state 
of plane polarization. However, the polarization gives more information than the curvature 
statistics. Indeed, zero curvature allows any structure that is constant along the magnetic 
field, in particular a set of magnetic filaments parallel to each other or a set of layered 
sheets, such that the magnetic field is constant on each sheet but its direction may change 
arbitrarily when passing from one layer to another. On the other hand, the wavenumber 
scalings obtained in this paper, (53) and (54), indicate that the magnetic field structures 
are layers in coordinate space and this rules out any filamentary structures. Further, our 
result about the plane polarization inhibit any "twists" of the magnetic field between layers, 
i.e. the magnetic field direction stays the same (or reverses) when passing from one layer to 
another. In fact, the presence of one neutral direction in the Lagrangian deformations tells 
us that these layers have a finite width in one direction and thus look like ribbons with the 
magnetic field directed along these ribbons. 

APPENDIX 

Our aim here is to derive a closed equation for the generating function Z (21), starting 
with equation (17). The last term in this equation is the easiest one 



where "D is a differential operator defined in (20). The correlators containing a factor of aij 
can be found using Gaussian integration by parts. In particular 



2«A;'((A|B(k)|' + Q;B'(k) + /3B'(k))£;) = -2«A;'PZ, 



(Al) 



{(TijE) = ) = —[X{{Tm,ijBjn + rm,ijBm)E) 



(A2) 



where we have used the definition (18). Here, F^.i^ is a response function 



F, 



5B, 



m 



(A3) 
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Differentiating (10) witfi respect to Uij and using tlie statistical wfiiteness of the strain tensor 
we get 



k,dj-^{l + kdi) 



(A4) 



In what follows we will make use of the isotropy of the turbulence, in particular, expressions 
of the type 



{B,B,E) 



2 " " |2 i7i\ / j: kikj 
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1 



\B\'E) S,, 



A;2 



(A5) 
(A6) 
(A7) 



Substituting (A4) into (A2) and using the above isotropy relations we have 



(A8) 



where "D is a differential operator defined in (20). This allows us to find the first term on 
the right hand-side of (17) 



kidj{aijE) 



^-j^k^Zkk + 3 (-2V + d^-l) kZk - 2VZ 
d d ^ ' 



(A9) 



Similarly, the other three terms on the right hand-side of (17) can be obtained via Gaussian 
integration by parts, and the use of the response function (A4) and isotropy condition. After 
some lengthy but straightforward algebra one obtains 



McTmliBmBl + BiBjn)E) = Afi 



(AlO) 



and 



2a{amiBmBiE) = 



+X{Zaf3 - Zxx) - ^kidiZx 



(All) 



-2/3Zap + 2/3Zxx - ^kAZa 
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The fourth term can be obtained from (All) via interchanging a with /3 and B with B 



2P{amiBmBiE) = pn 



-2aZap + 2aZxx - -^hdiZp 



(A12) 



Using the expressions (A9), (AlO), (All), (A12) and (Al), we obtain the final equation 



(A13) 



+ {2d - 6)VZ + A{1--)V'Z + 2{X' - 4al3){Z^f3 - Zxx) 



- 2Kk^VZ. 
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FIGURES 




FIG. 1. The magnetic field of 500 wavepackets in a numerical simulation with = 0.36 and 
K = 0. The figures show the particle's positions in A;-space at t = 6 and their corresponding real 
magnetic fields. The left and right figures correspond to different realisations of the strain matrix. 
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FIG. 2. The magnetic field of 500 wavcpackcts in a numerical simulation with 0, = 0.36 and 
K = 0.005. The figures show the particle's positions in k-sp&ce at t = 6 and their corresponding 
real magnetic fields for two different realisations of the strain matrix. The left figure corresponds 
to the same strain field as the left-hand graph of figure 1 and similarly for the right-hand graphs. 
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